Exact entropy of dimer coverings for a class of lattices in three or more dimensions 
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We construct a class of lattices in three and higher dimensions for which the number of dimer 
coverings can be determined exactly using elementary arguments. These lattices are a generalization 
of the two-dimensional kagome lattice, and the method also works for graphs without translational 
symmetry. The partition function for dimer coverings on these lattices can be determined also for 
a class of assignments of different activities to different edges. 
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The dimer model is a well-known problem in classical 
lattice statistics. The Ising model can be reformulated as 
a dimer model on a modified graph [l| . The dimer model 
is a simple, non-trivial but analytically tractable model 
which provides a starting point for the studies of systems 
of non-spherical molecules with hard-core interactions [2] ■ 
The latter can show a variety of geometrical phase tran- 
sitions e.g. in liquid crystals [SJ. Recently there has been 
a lot of interest in quantum dimer models, to describe 
the ground states of frustrated quantum antiferromag- 
nets, and some exotic phases expected in assemblies of 
hard core bosons on pyrochlore lattices [j| . 

Kasteleyn's pioneering solution of the problem for gen- 
eral two-dimensional planar lattices [5|, led to a lot of 
follow-up work on the two-dimensional problem [3|, but 
the problem in higher dimensions is much less studied. 
For non-zero fraction of sites not covered by dimers, it 
is known that there is nophase transition as a function 
of the density of dimers t(|| . Huse et al. have given non- 
rigorous arguments that the dimer-dimer orientational 
correlations have a power-law decay in all dimensions on 
bipartite lattices [7| . For graphs made of corner sharing 
triangles, where each vertex has a coordination number 
four, the entropy of dimer coverings per site can be deter- 
mined exactly in any dimension [8] . Huang et al studied 
the dimer problem in 3-dimensions, where the lattice is 
a stack of 2-d planes, but the problem is not fully 3- 
dimensional as the dimers were confined to lie in planes 
[9l| . Priezzhev calculated exactly the number of a subset 
of all dimer configurations on the cubic lattice ^Ifl] . For 
some anisotropic dimer models in three dimensions, the 
phase transition point can be determined exactly [ll|, 
but not the total entropy. 

Elser [l3| provided a simple argument to show that 
entropy per site of dimer covering of the kagome lattice 
(Fig. la) has a very simple form, and equals ilog2. 
This argument was extended to arbitrary graphs made 
of corner-sharing triangles by Misguich et al [gl- This 
problem was investigated again recently by Wang and 
Wu [I4I ■ The arguments of Misguich et al are valid only 
for graphs in which each vertex has coordination number 
4, and belongs to exactly two triangles. In this paper, we 
present an argument that is similar to the Elser's, but is 



valid for graphs in which the coordination number of ver- 
tices is not restricted to 4. Our arguments are applicable 
to a large class of lattices with arbitrary dimension, and 
also to graphs without translational invariance. Some 
examples of lattices where our method can be used are 
shown in Fig. 1. 

For simplicity of presentation, we first describe our ar- 
guments as applied to the kagome lattice (shown in Fig. 
la), and discuss the generalizations later. We select a 
subset of the sites of the lattice. We will call these the 
red sites, and denote the set by TZ. For the kagome lat- 
tice, we choose the red sites to be the sites of one of the 
three sublattices original lattice (Fig. la). For any dimer 
covering C of the lattice, to each red site s of the lattice, 
we assign a discrete variable ct^ . This variable takes only 
two possible values: CTs is -|-1 if the other vertex of the 
dimer covering s is above s, and —1 if it is below. 

We now arbitrarily choose a value of as for each site 
s in TZ. Let us denote this set by {cts}, and ask how 
many dimer configurations are consistent with {as}- We 
note that if CTs = +1, then the two bonds connecting 
the site s to sites below it would not be used in dimer 
coverings. Then we can safely delete these bonds, and 
consider dimer covering of the remaining graph. Simi- 
larly, if (Ts = —1, we remove two bonds connecting s to 
sites above. If we do this for all sites of the red sub- 
lattice, the full lattice breaks up into a set of mutually 
disconnected chains of the type shown in Fig. 2. 

It is easy to see that each chain can be covered by 
dimers in at most one way. For example, in Fig. 2, for 
the top chain, starting from left, we see that site 1 has to 
be matched with site 3, and then site 2 has to matched 
with 4. And so on. 

It is convenient to adopt the boundary conditions that 
a site on the right boundary can be either covered by 
a dimer, or left uncovered [14]. Then there is exactly 
one dimer covering consistent with any given set {cg}- 
The number of red sites is A^/3, and so the number of 
different possible choices of {as} is 2^^^. This is also the 
total number of dimer coverings. 

We define entropy per site C by the relation that num- 
ber of dimer coverings increases as exp(CiV), for a lattice 
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FIG. 1: Some examples of lattices for which the exact entropy is calculated in this paper, (a) The kagome lattice, (b) A 
three-dimensional lattice with corner-sharing triangles (c) The unit cell of the Na^Ir^Os lattice. There are 12 Iridium atoms 
per unit cell (shown in red or blue). Sodium and oxygen atoms are not shown. Atoms belonging to other unit cells are shown 
in dark red or black, (d) Unit cell of a lattice obtained by decoration procedure discussed in text. 



vith TV sites, for large TV. Then for the kagome lattice, 



C=-log2. 



(1) 



The above argument can be generalized to more gen- 
eral connected graphs made of corner-sharing triangles. 
Consider, for example, the lattices shown in Fig. lb and 
Ic. To construct the lattice of Fig. lb we start with a 
simple cubic lattice with only the vertical edges present. 
For each pair of neighbouring vertical chains we add a 
red site in between every fourth pair of opposing edges 
and join it to the ends of the edges. The opposing edges 
in between whom we add these additional sites are dis- 
placed with respect to each other so that no two triangles 
share an edge, as shown in Fig. lb. 

The unit cell of a hyperkagome lattice is shown in Fig. 
Ic. There are experimental materials like the sodium 
iridate (Na^Ir^Os) where this structure, known as the 



GGG structure, is realized in nature [l5|. To obtain this 
lattice we start from the pyrochlore lattice which con- 
sists of corner sharing tetragons, their centers forming 
a diamond lattice. The hyperkagome lattice in Fig. Ic 
is obtained by removing one site from each tetragon to 
leave a lattice of corner sharing triangles. 

For these cases, we choose the red sites as shown in 
Figs, lb and Ic, and as in the kagome case, to each red 
site assign a binary variable a, specifying whether the 
dimer covering that site belongs to one or the other of 
the two triangles meeting at that point. Then, as before, 
we ask how many different dimer coverings are consistent 
with an arbitrary choice oi {<Js}. 

If now, from each red site, we delete two edges that are 
not consistent with the given choice, the lattice breaks up 
into mutually disconnected chains of single edges and tri- 
angles similar to Fig. 2. These chains have a "backbone" 
of bonds in the z-direction in Fig. lb, and a zig-zag ar- 




FIG. 2: On deleting two edges from each red site, the kagome 
lattice breaks into disconnected chains. The values of the a 
variables at some of the red sites are indicated by numbers in 
parantheses. 



ray of roughly vertical lines in Fig. Ic (shown by black 
dashed lines). If we assume free boundary conditions on 
one (say the top) face, there is a unique dimer covering 
of the lattice consistent with any aribrary choice of {ps\- 
For large iV, the number of red sites is iV/3 (up to surface 
correction terms), and the number of dimer coverings is 
2^1'^ . Hence, for both the 3-dimensional lattices shown 
in Fig. lb and Ic, we get 
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For graphs in which each vertex is shared by two tri- 
angles, this result was obtained earlier by Misguich et al 
[31 . But our argument is extended easily to graphs where 
more than two triangles meet at a single vertex. Con- 
sider a graph formed by the vertices and edges of a col- 
lection of non-intersecting (possibly non-straight) lines, 
which may be imagined as embedded in a d-dimensional 
lattice (d > 2). We add a number of additional vertices 
not lying on any line, and call these additional vertices 
red vertices. The red vertices can be labelled by integers 
1, 2, . . . AV. For the jth red vertex, we choose a positive 
integer r^ > 1, choose Vj different edges from the original 
lines, and join both ends of each edge to the red vertex 
j by additional edges. Then the jth red vertex is the 
shared corner vetex of rj triangles, and its coordination 
number is 2rj . We ensure that the edges selected do not 
already have some vertex in common, and that no edge 
of the original lines is selected more than once, so that 
no edge belongs to more than one triangle. This defines 
the graph C for which we want to calculate the number 
of possible dimer covers. 

An example of such a graph is shown in Fig. Id. The 
lattice consists of vertices of a simple cubic lattice, la- 
belled by coordinates (mi,m2,?Ti3), where mi, 1712,1713 
are integers. We start with edges only along the z- 
direction. The lattice graph consists of disconnected ver- 



tical 1-dimensional chains. We add extra connections to 
this graph to make it 3-dimensional as follows: we pick 
an elementary cube of the lattice, add an extra vertex at 
the center of the cube, and connect it to two or more of 
the vertical edges of the cube by triangles (an example 
is shown in Fig. Id). We choose a unit cell of period 
pi,P2,P3 where pi,P2,P3 are integers > 1. (In the case 
shown in Fig. Id, Pi = P2 = P3 = 2.) We choose some 
extra cubes within the unit cell to connect this way, and 
then repeat this pattern to get a translationally invariant 
lattice. We ensure that these 'decorations' are selected 
such that the resulting graph is connected, and that no 
edge belongs to more than one triangle. 

In the example shown in Fig. Id, we chose two diago- 
nally opposite cubes within the unit cell, and add a site 
each at the center of each cube. The site with coordinates 
(1/2, 1/2, 1/2) is connected to all its four neighboring ver- 
tical edges, but we connect the site at (3/2, 3/2, 3/2) to 
only three of its neighboring vertical edges. 

Then, as before, at each red vertex j, we define a vari- 
able aj which takes rj possible values. There is a unique 
dimer cover of C consistent with a particular arbitrar- 
ily chosen set of values {tJs}- It follows that the total 
number of dimer covers of this graph Q is given by 



N,. 



(3) 



In the Fig. Id, we have two red vertices per unit cell, 
and ri = 4, and r2 = 3. There are 10 sites per cell. 
Hence the entropy per site for the lattice shown in Fig. 
Id is ^ log 12. 

It is easy to see that the arguments can be trivially ex- 
tended to other choices of decorations, and also to higher 
dimensional lattices. 

Our treatment can also be extended to the case where 
different edges are associated with different activities. 
This is illustrated most easily for the kagome lattice. Let 
us associate activity zi with a horizontal edge of the lat- 
tice (Fig. la). Let us assume that the activity for a 
non-horizontal edge is 2:2 if it belongs to an up-pointing 
triangle, and Z3 if it belongs to a down pointing trian- 
gle. Then, for a particular {cts} having n values -1-1, 
the weight of the unique dimer cover is easily seen to 
Z2 2:3 "-^1 ■ Summing over all choices of {cr}, we get 

total partition function ^{zi, Z2, Z3) = {z2 + Z3)^'^z-^ . 
This answer differs from that obtained by Wang and Wu 
[131, as our choice of edge weights differs from theirs. 

It is easy to see that this same method works for the 
higher dimensional graphs (e.g. for the lattice shown in 
Fig. Id), so long as (a) the two edges at a red vertex be- 
longing to the same triangle always have the same weight, 
and (b) all edges not having a red end vertex have the 
same weight. 

For example, consider the lattice shown in Fig. Id. 
In the unit cell shown in Fig. Id there are seven trian- 
gles. Let the non- vertical edges of the respective triangles 



have the activities zi, Z2, ■ ■ ■ z^ for the bonds in the 4 tri- 
angles meeting at (1/2, 1/2, 1/2), Zs, zg, zj for the bonds 
meeting at the site (3/2,3/2,3/2) and zq be the activity 
for the vertical edges. Then the free energy per site is 
(l/10)log[(zi + Z2 + Z3 + Z4){z5 + ze + zt)z^]. 

Note that the fact that a variables are completely un- 
correlated, suggests that orientational correlations be- 
tween dimers on these lattices may be short-ranged. In 
fact, for the kagome lattice, the dimer-dimer orienta- 
tional correlation function has been shown to be zero for 
all separations > 2 lattice spacings |12| . 

We should emphasize that though the {as} variables 
are local, independent variables, that uniquely specify 
the allowed dimer configurations, the procedure of finding 
the dimer configuration corresponding to a given {as} is 



nontrivial, and non-local. If we change as at only one red 
site, it changes the orientations of dimers very far along 
the two affected vertical chains. 
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